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Abstract. 

We introduce a generic Frohlich- Coulomb model of the oxides, which also includes 
infinite on-site (Hubbard) repulsion, and describe a simple analytical method of solving 
the multi-polaron problem in complex lattice structures. Two particular lattices, 
a zig-zag ladder and a perovskite layer, are studied. We find that depending on 
the relative strength of the Frohlich and Coulomb interactions these systems are 
either polaronic Fermi (or Luttinger)-liquids, bipolaronic superconductors, or charge 
segregated insulators. In the superconducting phase the carriers are super light mobile 
bipolarons. The model describes key features of the cuprates such as their T c values, 
the isotope effects, the normal state diamagnetism, pseudogap, and spectral functions 
measured in tunnelling and photoemission. We argue that a low Fermi energy and 
strong coupling of carriers with high-frequency phonons is the cause of high critical 
temperatures in novel superconductors. 
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There is overwhelming experimental |1|, |2|, [3|, [|, [5], and theoretical |7|, [S|, [9|, |K]] 

evidence for an exceptionally strong electron-phonon (el-ph) interaction in the cuprates, 
which competes with electron correlations. In recent years, several publications 
addressed the fundamental problem of competing el-ph and Coulomb interactions in 



the framework of the so-called Holstein-Hubbard model [D], [12|, 13, |14], I5(| , where both 
interactions are short-range (on-site). The model describes well many properties of 
the insulating state of the cuprates, including antiferromagnetism, lattice distortions, 
and phase segregation. However, it could hardly account for the high value of the 



superconducting critical temperature | 16fl . The mass of (bi)polaronic carriers in this 
model is very large in the relevant parameter region, and T c is suppressed below the 
kelvin scale. 

In choosing the correct interaction for HTSC, we take into account that most of 
the novel superconductors are doped insulators with highly polarizable ionic lattices. 
The low density of mobile carriers is unable to screen effectively the direct Coulomb 
electron-ion and electron-electron interactions. The layered structure of the cuprates 
reduces screening even further. Since the mobile carriers are confined to the copper- 
oxygen planes their interaction with out-of-plane ions, such as the apical oxygens, is 
particularly strong. A parameter-free estimate of the polaron binding energy E p in the 
cuprates puts it at about 0.5 eV or larger [17). This strong long-range (Frohlich) el-ph 
interaction necessarily leads to formation of small polarons. Such Frohlich small polarons 



were first considered in reference [[L8|]. Exact Monte-Carlo simulations of the single 
polaron problem [I9| showed that a long-range el-ph interaction effectively removes the 
difficulty with a large polaron (and bipolaron) mass in the Holstein-type el-ph models. 
Indeed, the polaron is heavy because it has to carry a lattice deformation with it which 
is the same deformation that forms the polaron itself. Therefore, there exists a generic 
relation between E p and the renormalization of its mass: m oc exp ( / yE p /u)), where uj is 
a characteristic phonon frequency and 7 ~ 1 is a numeric coefficient whose actual value 
depends on the radius of the interaction. For a short-range el-ph interaction (Holstein) 
the entire lattice deformation disappears and then forms in the new place when the 
polaron moves between the nearest lattices sites. Therefore, 7 = 1 and the polaron is 
very heavy for the characteristic cuprates values E p ~ 0.5 eV and uj ~ 0.05 eV. In the 
case of a long-range interaction, only a fraction of the total deformation changes every 



time the polaron moves and 7 could be as small as 0.25 pq . Clearly, this results in 



a dramatic mas reduction since 7 enters the exponent. Thus the effective mass could 
be < 10 m e where a naive Holstein-like estimate would yield a huge mass ~ 10,000m e . 
The above qualitative reasoning was fully confirmed by analytical [0 and numerical 
(approximation- free Monte Carlo) |19j studies of the double-chain and double-plane 
models with long-range el-ph interactions. Later the single-polaron and bipolaron cases 
of the chain model were analyzed in more detail in references [^(J and [^1|], respectively. 
These studies confirmed a much lower mass of both polaron and bipolaron in comparison 
to the Holstein-Hubbard limit. 

Here we argue that a consistent theory of HTSC must include both the long- 
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range Coulomb repulsion between the carriers and the strong long-range electron- 
phonon interaction. We propose an analytically solvable multi-polaron model of high- 
temperature superconductivity that includes these realistic long-range interactions. 
From theoretical standpoint, the long-range Coulomb repulsion is critical in ensuring 
that the carriers would not form large clusters. Indeed, in order to form stable pairs 
(bipolarons) the el-ph interaction has to be strong enough to overcome the Coulomb 
repulsion at short distances. Since the el-ph interaction is long-range, there is a potential 
possibility for clustering. We shall demonstrate that the inclusion of the Coulomb 
repulsion V c makes the clusters unstable. More precisely, there is a certain window of 
V c /Ep inside which the clusters are unstable but bipolarons nonetheless form. In this 
parameter window the bipolarons are light and the system is a superconductor with a 
high critical temperature. The bipolarons repel each other and propagate in a band of 
about the same bandwidth as the single-polaron bandwidth, in sharp contrast with all 
bipolaronic models considered previously. At a weaker Coulomb interaction the system 
is a charge segregated insulator. At a stronger Coulomb repulsion the system is a polaron 
Fermi (or Luttinger) liquid. In the superconducting phase but close to the clustering 
boundary, dynamical formation of short lived clusters or stripes could be expected. 

Our generic Frohlich-Coulomb model explicitly includes the electron kinetic energy, 
the infinite-range Coulomb and electron-phonon interactions as well as the lattice energy. 
The implicitly present infinite Hubbard U prohibits double occupancy and removes the 
need to distinguish the fermionic spin. Introducing spinless fermion operators c n and 
phonon operators d ma , the model Hamiltonian is written as 

H = - T ( n _ n ') C n C n' + V d n ~ nO^Pni/Cn' 

n^n' n^n' 

- co^2g a (m - n)(e 

ma ' U m _ n )c n C n ((i mQ , -|- d ma 
nm ma ^ ' 

We note that the el-ph term is written in real rather than momentum space. This is 
more convenient in working with complex lattices. Here e ma is the polarization vector 
of ath vibration coordinate at site m, u m _ n = (m — n)/|m — n| is the unit vector in 
the direction from electron n to the ion m, and g a (m — n) is the dimensionless el-ph 
coupling function. [g a (m — n) is proportional to the force acting between m and n.] 
We assume that all the phonon modes are dispersionless with frequency u and that the 
electrons do not interact with displacements of their own atoms, g a (0) = 0. We also use 
h = 1 throughout the paper. 

In general, the many-body model (1) is of considerable complexity. However, we 
are interested in the limit of strong el-ph interaction. In this case, the kinetic energy is 
a perturbation and the model can be grossly simplified in a two-step procedure. On the 
first step, the Lang-Firsov canonical transformation is performed which diagonalizes 
the last three terms in equation (1). Introducing S = Emn«fi l a( m ~ n )( e m« • 

Um- n)c n C n 

(d-lna ~ d ma ) one obtains the transformed Hamiltonian without an explicit 
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el-ph term 

H = e~ s He s = - ^ 0w4 c n' + ^Y. (^Lo^mo + ~ 

n^n' ma V 

+ v(n - n^C^CnC^Cn' - E p Y 4i C n- (2) 

n^n' n 

The last term describes the energy which polarons gain due to el-ph interaction. E p is 
the familiar polaron (Franc-Condon) shift 

E p = w^^(m - n)(e mQ • u m _ n ) 2 , (3) 

mo 

which we assume to be independent of n. E p is a natural measure of the strength of the 
el-ph interaction. The third term in Eq.(2) is the polaron-polaron interaction: 

v(n - n') = V c (n - n') - V pa (n - n'), (4) 
Vp a (n - n') = 2u Y 9a(m - n)# Q (m - n') x 

mo 

(^ma ' "m- n)(^mo ' U m _ n '), (5) 



where V pa is the inter-polaron attraction due to joint interaction with the same vibrating 
atoms. Finally, the first term in equation (2) contains the transformed hopping operator 

Onn': 



Onn' 



T(n — n') exp 



Y [g a {m - n)(e ma • u m _ n ) 

_mo 

- g a (m-n')(e ma - u m ^ n/ )}(dl na - d ma ) . (6) 

At large E p /T(n — n') this term is a perturbation. In the first order of the strong 
coupling perturbation theory @, cr nn / should be averaged over phonons because there 
is no coupling between polarons and phonons in the unperturbed Hamiltonian [the last 
three terms in equation (2)]. For temperatures lower than u, the result is 

f(n - n') = (a nn ,) ph = T(n - n') exp[-G 2 (n - n')], (7) 
G 2 (n - n') = Y 9a (m - n) (e mQ ■ u m _ n ) x 

ma 

[g a (m - n)(e mQ • u m _ n ) - g a (m - n')(e ma ■ u m _ n /)] . (8) 

By comparing equations (3) and (5) with equation (8), the mass renormalization 
exponents can be expressed via E p and V pSu as follows 

G 2 (n-n') = ^(i? p -^ pa (n-n')). (9) 

This results in a renormalized hopping term that represents the small parameter of 
a strong coupling perturbation theory |7[]. The above technical transformation is 
simple and has been described elsewhere together with a detailed description of 
the perturbation procedure. The resulting model is purely polaronic in which phonons 
are "integrated out" 

H p = H + f/pert, (10) 
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H = -E p C n C n + Y V ( n ~ n ') C i C n C i> C n>, (H) 
n n ^ n ' 

-ffpert = - Y t ( n ~ n ') C n C n'- (12) 
n^n' 

When Vp a exceeds V c the full interaction becomes negative and polarons form pairs. 
We emphasize that while the above formalism fails in some regimes of the short-range 
el-ph models (for instance, the adiabatic limit of the Holstein model), it is surprisingly 
accurate for long-range el-ph interactions, as was demonstrated by comparing the 
analytical results with exact quantum Monte Carlo data (T^j. It makes theoretical 
analysis of even complex interactions and lattice geometries simple and instructive. But 
before proceeding to analyzing concrete lattices, let us elaborate on the physics behind 
the lattice sums in equations (3) and (5). 

When a carrier (electron or hole) acts on an ion with a force f , it displaces the ion 
by some vector x = f/s. Here s is the ion's force constant. The total energy of the 
carrier-ion pair is — f 2 /(2s). This is precisely the summand in equation (3) expressed 
via dimensionless coupling constants. Now consider two carriers interacting with the 
same ion, see figure 1(a). The ion displacement is x = (fj + f 2 )/s and the energy is 
— if /(2s) — ff/(2s) — (fi ■ f2)/s. The last term here should be interpreted as an ion- 
mediated interaction between the two carriers. It depends on the scalar product of fi 
and f 2 and consequently on the relative positions of the carriers with respect to the 
ion. If the ion is an isotropic harmonic oscillator, as we assume in this paper, then 
the following simple rule applies. If the angle </> between fi and f 2 is less than 71/ 2 
then the polaron-polaron interaction is attractive, otherwise it is repulsive, see figure 
1(b). The overall sign and magnitude of the interaction is given by the lattice sum in 
equation (5), evaluation of which is elementary. Notice also that according to equation 
(9), an attractive interaction reduces the polaron mass (and consequently bipolaron 
mass), while repulsive interaction enhances the mass. Thus in our model the long- 
range character of the el-ph interaction serves the double purpose. Firstly, it generates 
additional inter-polaron attraction because the distant ions have small angle (p. This 
additional attraction helps overcome the direct Coulomb repulsion between the polarons. 
Secondly, the Frohlich interaction makes the bipolarons light leading to a high critical 
temperature. 

The many-particle ground state of H depends on the sign of the polaron-polaron 
interaction, the carrier density, and the lattice geometry. First we consider the zig- 
zag ladder, figure 2(a), assuming that all sites are isotropic two-dimensional harmonic 
oscillators. For simplicity, we also adopt the nearest-neighbour approximation for both 
interactions, g a Q) = 9, V c (n) = V c , and for the hopping integrals, T(m) = T/vjv > for 
I = n = m = a, and zero otherwise. Hereafter we set the lattice period a — 1. There 
are four nearest neigbours in the ladder, z = 4. The one-particle polaronic Hamiltonian 
takes the form 

Hp = - Y ( E p[4i c n + PlPn] + t'[c j n+1 c n + pl +1 p n + h.c] + t[p ] n c n + p\_ x c n + h.c.]) , (13) 

n 
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where c n and p n are polaronic operators on the lower and upper side of the ladder, 
respectively, see figure 2(b). Applying the general formulas (3), (5), and (9), we 
obtain E p = Ag 2 u, t' = T^n exp[— 7E p /(Su)], and t = T^n exp[— 3E p /(Au)]. Fourier 
transformation yields the one-particle spectrum 

E 1 {k) = -E p - 2t'cos(fc) ±icos(fc/2). (14) 

Two overlapping polaronic bands have a combined width of W = At' + It. The lower 
band has the bandwidth W and the effective mass m\ = 2 /(At' + 1) near the bottom, 
while the upper band has the bandwidth At' — 2t and a heavier mass m* = 2/ (At' — t). 

Let us now place two polarons on the ladder. The nearest neighbour interaction, 
equation (4), is found as v — V c — E p /2 if two polarons are on the different sides of 
the ladder, and v = V c — E p /A if both polarons are on the same side. The attractive 
interaction is provided via the displacement of the lattice sites which are the common 
nearest neighbours to both polarons, under the condition that the angle between the 
directions pointing from those sites to two polarons is less than tt/2. If > 7r/2 the 
effective interaction of two polarons is repulsive. There are two such nearest neighbours 
for the intersite bipolaron of the type A or B, figure 2(c), but there is only one common 
nearest neighbour for the bipolaron C, figure 2(d). When V c > E p /2, there are no bound 
states and the multi-polaron system is a ID Luttinger liquid. However, when V c < E p /2 
and consequently v < 0, the two polarons are bound into an intersite bipolaron of the 
type A or B. 

It is quite remarkable that the bipolaron tunnelling appears already in the first order 
in polaron hopping H pcrt as was anticipated in ||16|| . This case is different from both the 
on-site bipolaron discussed a long time ago ||23|| , and from the inter-site chain bipolaron 
discussed recently ^TJ, where the bipolaron tunnelling was of the second order in t. 
Indeed, in the first order in H pert one should consider only the lowest energy degenerate 
configurations A and B and discard the processes that involve all other configurations. 
The result of such a projection is a bipolaronic Hamiltonian 

H b = (V c - ^-E p ) £[44. + B\B n \ - t'Y\BlA n + B\_ X K + h.c], (15) 



2 

where A n = c n p n and B n = p n c n+ \. Fourier transformation yields the bipolaron energy 
spectrum: 

E 2 (k) = V c -^E p ±2t' cos(k/2). (16) 

There are two bipolaron bands with a combined width of At'. The bipolaron binding 
energy is 

A = 2E 1 (0) - E 2 (0) = ^ - V c -2t-At'. (17) 

The bipolaron mass near the bottom of the lowest band is m** = 2/t'. Neglecting t and 
t' relative to E p and V c we arrive at the following conclusion. When V c < E p /2, two 
polarons form a bipolaron with effective mass m** « (4 + exp^)m^. The numerical 
coefficient | ensures that m** remains of the order of m* even at large E p . 
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In models with strong intersite attraction there is a possibility of clasterization. In 
a way similar to the two-particle case above, the lowest energy of n polarons placed on 
the nearest neighbours of the ladder is found as E n = (2n — 3)V C — ^f^-E p , for any n > 3. 
There are no resonating states for n-polaron nearest neighbour configuration if n > 3. 
Therefore there is no first-order kinetic energy contribution to their energy. E n should 
be compared with the energy E\ + (n — l)E 2 /2 of far separated (n — l)/2 bipolarons 
and a single polaron for odd n > 3, or with the energy of far separated n bipolarons for 
even n > 4. "Odd" clusters are stable at V c < 6n ™ 1Q E p , and "even" cluster are stable at 
V c < ^zj^Ep. Here we have neglected the kinetic energy of polarons and bipolarons. As 
a result we find that bipolarons repel each other and single polarons at V c > ^E p . If V c is 
less than ^E p then immobile bound clusters of three and more polarons could form. We 
would like to stress that at distances much larger than the lattice constant the polaron- 



polaron interaction is always repulsive [16|, and the formation of infinite clusters, stripes 



or strings is impossible p4| . Combining the condition of bipolaron formation and that 



of the instability of larger clusters we obtain a window of parameters 

h p <V C < l -E P) (18) 

within which the ladder is a bipolaronic conductor. Outside this window the ladder is 
either a charge segregated insulator (small V c ) or the one-dimensional (ID) Luttinger 
liquid (large V c ). 

Our consideration is directly related to doped cuprates. Here we consider a two 
dimensional lattice of ideal octahedra that can be regarded as a simplified model of 
the copper-oxygen perovskite layer, figure 3. The lattice period is a = 1 and the 
distance between the apical sites and the central plane is h = a/2 = 0.5. All in- 
plane atoms, both copper and oxygen, are static but apical oxygens are independent 
three-dimensional isotropic harmonic oscillators. Because of poor screening the hole- 
apical interaction is purely Coulombic, g a (m — n) = n a /\m — n| 2 , a = x,y,z. To 
account for the experimental fact that z-polarized phonons couple to the holes stronger 
than the others 0, we choose /\/2- The direct hole-hole repulsion is 

V c (n — n') = so that the repulsion between two holes in the NN configuration is 

V c . We also include the bare NN hopping T/vjv; the next nearest neighbor (NNN) hopping 
across copper T NNN , and the NNN hopping between the pyramids T' NNN . According 
to equation (3), the polaron shift is given by the lattice sum (after summation over 
polarizations): 

E p = 2k*o;£ ( 1 + j h " ) = 31.154a,, (19) 
^\Jm — n| 4 |m — n) / 

where the factor 2 accounts for the two layers of apical sites. [For reference, Cartesian 
coordinates are n = (n x + 1/2, n y + 1/2,0), m = (m x ,m y , h); n x ,n y ,m x ,m y being 
integers.] The polaron-polaron attraction is 

. . 9 ^ h 2 + (m — n') • (m — n) . . 

V pa n - n' = 4w4 £ 20 

p „ m-n /3 m-n 3 
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Performing the lattice summations for the NN, NNN, and NNN' configurations one 
finds Vp a = 1.23 E p , 0.80 E p , and 0.82 E p , respectively. Substituting these results in 
equation (4) and (9) we obtain the full inter-polaron interaction: vnn = V c — 1.23 E p , 
vnnn = tt| — 0.80 -Ep, v' NNN = — 0.82 E p , and the mass renormalization exponents: 
G 2 NN = 0.38(E p /co), G 2 NNN = 0.Q0(E p /uj) and G% NN = 0.59(E p /lu). 

Let us now discuss different regimes of the model. At V c > 1.23 E p , no bipolarons 
are formed and the systems is a polaronic Fermi liquid. The polarons tunnel in 
the square lattice with NN hopping t = T NN exp(— 0.38E p /u) and NNN hopping 
t' = T NNN exp(— 0.Q0E p /uj). [Since G 2 NNN « G% NN one can neglect the difference 
between NNN hoppings within and between the octahedra.] The single polaron spectrum 
is therefore 

#i(k) = -E p - 2t'[cosk x + cosy - At cos(k x /2) cos(k y /2). (21) 

The polaron mass is m* = l/(t + 2t'). Since in general t > t', the mass is mostly 
determined by the NN hopping amplitude t. If V c < 1.23 E p then intersite NN bipolarons 
form. The bipolarons tunnel in the plane via four resonating (degenerate) configurations 
A, B, C, and D, see figure 4. In the first order in H pcrt one should retain only these 
lowest energy configurations and discard all the processes that involve configurations 
with higher energies. The result of such a projection is a bipolaron Hamiltonian 

H b = (V c - 3.23 E p ) Y\A A \ + B\B X + C}d + D\Di] 

i 

- t'Y\ABi + Bid + C\D\ + D\A, + h.c] 

t'Y\A\_^B x + B} +y G + Cf +X A + D\_ y Ai + h.c], (22) 

n 

where 1 numbers octahedra rather than individual sites, x = (1,0), and y = (0,1). 
A Fourier transformation and diagonalization of a 4 x 4 matrix yields the bipolaron 
spectrum: 

E 2 (k) =V C - 3.23E P ± 2t'[cos(Jfe a ./2) ± cos(ib v /2)]. (23) 

There are four bipolaronic subbands combined in the band of the width 8t'. The 
effective mass of the lowest band is m** = 2/t'. The bipolaron binding energy is 
A = 1.23E P — V c — A(2t + t'). As in the ladder, the bipolaron moves already in the 
first order in polaron hopping. This remarkable property is entirely due to the strong 
on-site repulsion and long-range electron-phonon interaction that leads to a non-trivial 
connectivity of the lattice. This situation is unlike all other models studied previously. 
[Usually, the bipolaron moves only in the second order in polaron hopping and therefore 
is very heavy] In our model, this fact combines with a weak renormalization of t! 
yielding a superlight bipolaron with mass m** oc exp(0.60 E p juS). We recall that in the 
Holstein model m** oc exp(2E p /uj). Thus the mass of the Frohlich bipolaron scales 
approximately as a cubic root of that of the Holstein one. 

At even stronger el-ph interaction, V c < l.\QE P) NNN bipolarons become stable. 
More importantly, holes can now form 3- and 4-particle clusters. The dominance of 
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the potential energy over kinetic in Hamiltonian (10) enables us to readily investigate 
these many-polaron cases. Three holes placed within one oxygen square have four 
degenerate states with energy 2{V C — 1.23-Ep) + ^ — 0.802^. The first-order polaron 
hopping processes mix the states resulting in a ground state linear combination with 
energy E 3 = 2.71V C — 3.2QE P — V 'At 2 + t' 2 . It is essential that between the squares such 
triads could move only in higher orders in polaron hopping. In the first order, they are 
immobile. A cluster of four holes has only one state within a square of oxygen atoms. 
Its energy is £ 4 = 4(V C - 1.23E P ) + 2(^ - 0.80E P ) = 5A1V C - Q.52E p . This cluster, as 
well as all the bigger ones, are also immobile in the first order of polaron hopping. We 
conclude that at V c < 1.16 E p the system quickly becomes a charge segregated insulator. 

The fact that within the window 1.1QE P < V c < 1.23E P there are no three or higher 
polaron bound states, means that bipolarons repel each other. The system is effectively 
the charged Bose-gas which is a well known superconductor 0. The superconductivity 
window that we have found, is quite narrow (see figure 5). This indicates that the 
superconducting state in such systems is a subtle phenomenon which requires a fine 
balance between electronic and ionic interactions. Too strong el-ph interaction leads to 
clustering, while too weak interaction cannot bind the carriers and the superconductivity 
is at best of BCS type. These considerations may provide additional insight into the 
uniqueness of one particular structure, the copper-oxygen perovskite layer, for HTSC. 
It also follows from our model that superconductivity should be very sensitive to any 
external factor that affects the balance between V c and E p . For instance, pressure 
changes the octahedra geometry and hence E p and V v& . Chemical doping enhances 
internal screening and consequently reduces E p . 

We now assume that the superconductivity condition is satisfied and show that our 
Frohlich- Coulomb model possesses many key properties of the underdoped cuprates. 
The bipolaron binding energy A should manifest itself as a normal state pseudogap 
with size of approximately half of A [[?]]. Such a pseudogap was indeed observed in many 
cuprates. In contrast with the BCS superconductor, the symmetry of the pseudogap 
might differ from the symmetry of the superconducting order parameter, which depends 
on the bipolaronic band dispersion. The symmetry of the order parameter was found 
to be <2-wave f25j, while the former is an anisotropic s-wave, in accordance with many 
experimental observations. There should be a strong isotope effect on the (bi)polaron 
mass because t, t' oc exp(— const \fM). Therefore the replacement of O 16 by O 18 increases 



the carrier mass [ 2"6f . Such an effect was observed in the London penetration depth of 
the isotope-substituted samples [l|]. The mass isotope exponent, a m = d In m**/d In M, 
was found to be as large as a m = 0.8 in Lai. 8 9 5 Sr .i05CuO4. Our theoretical exponent 
is a m = 0.3E p /u, so that the bipolaron mass enhancement factor is exp(0.6E p /uj) ~ 5 
in this material. With the bare hopping integral T^nn = 0.2 eV we obtain the in- 
plane bipolaron mass m** ~ 10 m e . Calculated with this value the in-plane London 
penetration depth, X a b = [m**/87rne 2 ] 1 / 2 ~ 316 nm (n the hole density) agrees well with 
the measured one \ a b — 320 nm. Taking into account the c-axis tunnelling of bipolarons, 
the critical temperature of their Bose-Einstein condensation can be expressed in terms 



Frohlich- Coulomb model of high-temperature superconductivity 



10 



of the experimentally measured in-plane and c-axis penetration depths, and the in-plane 
Hall constant Rh as T c « 1.64(e.Rtf/A^A;?) 1//3 . Here T c , eRn, and A are measured in K, 
cm 3 and cm, respectively ||27|| . Using the experimental X a b = 320 nm, A c = 4160 nm, 
and Rh = 4 x 10~ 3 cm 3 /C (just above T c ), one obtains T c = 31 K in striking agreement 
with the experimental value T c = 30 K. The recent observation of the normal state 
diamagnetism in La 2 _a ; Sr a .Cu04 [^8] also fits well the prediction of the bipolaron theory 
29fl . Many other features of the bipolaronic (super) conductor, e.g., the unusual upper 



critical field, electronic specific heat, optical, ARPES and tunnelling spectra match those 

in the cuprates (for a recent review, see reference |30| ) . 

Finally, we show that the Fermi energy in all novel superconductors is surprisingly 

low, of the order or even smaller than the most essential optical phonon energy. The 

band structure of the cuprates is quasi-two-dimensional with a few degenerate hole 

pockets. Applying the parabolic approximation for the band dispersion one obtains the 

renormalized (polaronic) Fermi energy as 

7m,- d , , 

Ef = — b 24 

where d is the interplane distance, and rii, m* are the density of holes and their effective 
mass in each of the hole subbands i renormalized by the electron-phonon (and electron- 
electron) interactions. One can express the renormalized band-structure parameters 
through the in-plane London penetration depth at T = 0, measured experimentally: 

A H i m i 

As a result, one obtains the parameter-free expression for the Fermi energy as 
4#e 2 A^ 

where g is the degeneracy of the spectrum. The degeneracy g in the cuprates may depend 
on doping. In underdoped cuprates one expects 4 hole pockets inside the Brillouin 
zone (BZ) due to the Mott-Hubbard gap. If the hole band minima are shifted with 
doping to BZ boundaries, the spectrum will be two-fold degenerate, so that g > 2 in 
cuprates. Because equation (|26|) does not contain any other band-structure parameters, 
the estimate of Ep using this equation does not depend very much on the parabolic 
approximation for the band dispersion. Generally, the ratios n/m in equations (|2~4]) and 
(f25|) are not necessary the same. The "superfluid" density in equation ( p5|) might be 
different from the total density of delocalized carriers in equation (pip. However, in 



a translational invariant system they must be the same This is true even in the 
extreme case of a pure two-dimensional superfluid, where quantum fluctuations might be 
important. One can obtain a reduced value of the zero temperature superfluid density 
only in the dirty limit I <C £(0) where £(0) is the zero-temperature coherence length. 
The latter was measured directly in cuprates as the size of the vortex core. It is about 
10 A or even less. On the contrary, the mean free path was found surprisingly large at 
low temperatures, I ~ 100-1000 A. Hence, the cuprates are in the clean limit, / ^> £(0), 
so that the parameter-free expression for Ep, equation (pB]), is perfectly applicable. 
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Equation (26) yields < 100 meV for the cuprates, especially if the degeneracy 
g > 2 is taken into account. A few examples are La1.85Sro.15Ci1.O4 (T c = 37 K, X H = 240 
nm 0, d = 0.66 nm) with gE F = 77 meV, YBa 2 Cu 3 6 .92 (T c = 91.5 K , \ H = 186 
nm f2"? |, d = 0.43 nm) with gEp = 84 meV. That should be compared with the 
characteristic phonon frequency, which is estimated as the plasma frequency of oxygen 
ions, u = (AitZ 2 e 2 N/M) 1 / 2 . One obtains u = 84 meV with Z = 2, N = 6/V cell , M = 16 
a.u. for YBa2Cu30g. Here V ce u is the volume of the chemical unit cell. The low Fermi 
energy Ep < ui is a serious problem for the Migdal-Eliashberg approach. The non- 
crossing diagrams cannot be treated as vertex corrections because u/Ep > 1, since they 
are comparable to the standard terms. On the contrary, the estimate of Ep supports 
further the nonadiabatic (bi)polarons as the (super) carriers in high-T c superconductors. 

In conclusion, we have introduced a realistic multi-polaron model of high- 
temperature superconductivity with the strong Frohlich and Coulomb long-range 
interactions. We have described a simple procedure of calculating polaron and bipolaron 
masses, identified and quantitatively analyzed a new resonance mechanism of bipolaron 
mass reduction, and found the conditions for clustering of holes and the window for 
their high-T c superconductivity. The model possesses a rich phase diagram in the 
coordinates of the inter-site Coulomb repulsion V c and the polaronic (Franck-Condon) 
level shift E p , see figure 5. The ground state is a polaronic Fermi (or Luttinger) 
liquid, for the strong Coulomb repulsion, bipolaronic high-temperature superconductor 
for the intermediate Coulomb repulsion, and the charge-segregated insulator for the 
weak repulsion. Remarkably, the inter-site bipolarons in the superconducting phase are 
"superlight" , propagating coherently with about the same mass as single polarons. In 
our model the bipolarons tunnel already in the first order in polaron tunnelling which 
results in the bipolaron mass scaling linearly with the polaron hopping integral. Many 
properties of the model in the superconducting phase match those of the cuprates. We 
argue that a surprisingly low Fermi energy and the strong unscreened coupling of carriers 
with high frequency optical phonons is the origin of high temperature superconductivity. 

This work has been supported by EPSRC UK (grant R46977), by the Leverhulme 
Trust (grant F/00261/H), and by DARPA. 
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Figure 1. The mechanism of the polaron-polaron interaction, (a) Together, the two 
polarons (solid circles) deform the lattice more effectively than separately. An effective 
attraction occurs when the angle 4> between xi and X2 is less than ir/2. (b) A mixed 
situation. Atom 1 results in repulsion between two polarons while atom 2 results in 
attraction. 
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Figure 2. Onc-dimcnsional zig-zag ladder, (a) Initial ladder with the bare hopping 
amplitude T. (b) Two types of polarons with their respective deformations, (c) 
The two degenerate bipolaron configurations A and B. (d) A different bipolaron 
configuration C which energy is higher than that of A and B. 




Figure 3. A fragment of the perovskite layer. 




Figure 4. Four degenerate bipolaron configurations A, B, C, and D. Some single- 
polaron hoppings are indicated by arrows. 



Frohlich- Coulomb model of high-temperature superconductivity 17 




insulator 



Polaron shift E p 



Figure 5. Phase diagram of the Frohlich- Coulomb model. The model is a 
polaronic Fermi liquid for the strong Coulomb repulsion, bipolaronic high-temperature 
superconductor (HTSC) for the intermediate Coulomb repulsion, and the charge- 
segregated insulator for the weak repulsion. 



